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Kaon production in pion-nucleon collisions in nuclear mat- 
ter is studied in the resonance model. To evaluate the in- 
medium modification of the reaction amplitude as a function 
of the baryonic density we introduce relativistic, mean-field 
potentials for the initial, final and intermediate mesonic and 
baryonic states. These vector and scalar potentials were cal- 
culated using the quark-meson coupling (QMC) model. The 
in-medium kaon production cross sections in pion-nucleon in- 
teractions for reaction channels with A and S hyperons in the 
final state were calculated at the baryonic densities appro- 
priate to relativistic heavy ion collisions. Contrary to earlier 
work which has not allowed for the change of the cross sec- 
tion in medium, we find that the data for kaon production are 
consistent with a repulsive Jf^-nucleus potential. 



I. INTRODUCTION 

The properties of kaons in nuclear matter have re- 
cently attracted enormous interest because of their ca- 
pacity to signal chiral symmetry restoration or give infor- 
mation on the possibility of kaon condensation in neutron 
stars [0-D . Studies with a variety of models [||-|| indicate 
that the antikaon potential is attractive while the kaons 
feel a repulsive potential in nuclear matter. The results 
from kaonic atoms , as well as an analysis [TTl-[r^ ] 
of the K~ production from heavy ion collisions [[l6pp0f, 
are in reasonable agreement with the former expectation 
for antikaons. However, the analysis of available data 
on K'^ production from heavy ion collisions at SIS en- 
ergies ll^-^ contradicts the predictions that the kaon 
potential is repulsive. The comparison between the heavy 
ion calculations and the data p"3|- p^ , ^ , p2[ indicates that 
the i4r+-meson spectra are best described by neglecting 
any in-medium modification of the kaon properties. Fur- 
thermore, the introduction of even a weakly repulsive 
if^-nucleus potential results in a substantial underesti- 
mate of the experimental data on kaon production. 

Since in heavy ion collisions at SIS energies 16-20|] 
the A'+-mesons are predominantly produced by sec- 
ondary pions, we investigate the kaon production reac- 
tions, TT^N^Y+K (y = A, E hyperons), in nuclear mat- 
ter. To be specific, we combine earlier studies of kaon 
production in free space with a very successful, relativis- 
tic mean field description of nuclear systems (QMC). All 
parameters are fixed by earlier studies and the effects 



of the medium on the reaction cross sections are calcu- 
lated for the first time. The result is impressive in that 
the medium effects explain the nuclear production data, 
without any adjustment of the parameters determined 
elsewhere, including the standard repulsive kaon-nucleus 
interaction. 

Our paper is organized as follows. In Sec. II we in- 
troduce the vector and scalar potentials for mesons and 
baryons involved in the calculations of the tt+N^Y-^-K 
amplitudes. We explain in Sec. Ill the resonance 
model which is used to calculate the cross sections, 
TT+N^Y+K . The strangeness production threshold in 
nuclear matter and its dependence on the baryon density 
is discussed in Sec. IV. The cross sections of kaon produc- 
tion in TT+N collisions in vacuum and in nuclear matter 
at different densities are then evaluated and shown in 
Sec. V for the Tr~ +p^A+K^ reaction and in Sec. VI for 
the TT+N^Y^+K reactions. The impact of our results on 
heavy ion collisions is discussed in Sec. VII. Finally, the 
summary and conclusions are given in Sec. VIII. 



II. MEAN-FIELD POTENTIALS FOR MESONS 
AND BARYONS 



In the present study, we use the quark-meson cou- 
pling (QMC) model ||2^, which has been successfully 
applied not only to the problems of conventional nu- 
clear physics [ p4|j25| but also to the studies of meson 
and hyperon properties in a nuclear medium |]7|,p6|-^. 
A detailed description of the Lagrangian density and the 
mean- field equations are given in Refs. p|p4[ ^j3l| ]. The 
Dirac equations for the quarks and antiquarks in the 
hadron bags {q = u,u,d or d, hereafter), neglecting the 
Coulomb force, are given by: 



V'd(a;) 



[i^ ■ dx - TOs] ipsix) (or ipsix)) 0. 



(1) 



(2) 
(3) 



The mean-field potentials for a bag in nuclear matter 
are defined by V.^^g^^a, VJ= g^co and Vp^=g'/,b, with g$, 
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and the corresponding quark-meson coupling con- 
stants. 

The normaUzed, static sohition for the ground state 
quarks or antiquarks in the hadron, h, may be written 



as 



(4) 



where / labels the quark flavours, and Nf and tpfix) 
are the normalization factor and corresponding spin and 
spatial part of the wave function. The bag radius in 
medium, i?^, which generally depends on the hadron 
species to which the quarks and antiquarks belong, will 
be determined through the stability condition for the (in- 
medium) mass of the meson against the variation of the 
bag radius @j2j-§^,0 (see also Eq. (||)). The eigenen- 
ergies, £/, in Eq. Ij) in units of are given by 



where ^l* = 



a;2 + (i?*m;)2, with 
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(5) 

(6) 
(7) 

and 
and 

X g , are determined by the usual, linear boundary con- 
dition |2^,^. Note that the lowest eigenenergy for the 
Dirac equation (Hamiltonian) for the quark, which is pos- 
itive, can be regarded as the analog of a constituent quark 
mass. 

The hadron masses in symmetric nuclear matter rele- 
vant for the present study are calculated by: 



ils — \/x^+iRfjns)^. The bag eigenfrequencies. 



{Uq + nq)W + {Us + ns)Vls - Zh 



Rl 



dmt 



dRh 



-TTRfB, (8) 



(9) 



In Eq. I 



(riq) and Us (rij) are the lowest mode 



light quark (antiquark) and strange (antistrange) quark 
numbers in the hadron, h, respectively, and the Zh 
parametrize the sum of the ccnter-of-mass and gluon fluc- 
tuation effects, and are assumed to be independent of 
density. The parameters are determined in free space to 
reproduce their physical masses. 

In this study we chose the values niq—b MeV and 
■ms—2bQ MeV for the current quark masses, and Rn=0.8 
fm for the bag radius of the nucleon in free space. Other 
input parameters and some of the quantities calculated 
are given in Refs. p^-^. We stress that while the model 
has a number of parameters, only three of them, g^, g^ 
and g'^, are adjusted to fit nuclear data - namely the sat- 
uration energy and density of symmetric nuclear matter 
and the bulk symmetry energy. None of the results for 
nuclear properties depend strongly on the choice of the 



other parameters - for example, the relatively weak de- 
pendence of the final results for the properties of finite 
nuclei, on the chosen values of the current quark mass 
and bag radius, is shown explicitly in Refs. p^ , ^ . Ex- 
actly the same coupling constants, g^, g^ and g^, are 
used for the light quarks in the mesons and hyperons as 
in the nucleon. However, in studies of the kaon system, 
we found that it was phenomenologically necessary to 
increase the strength of the vector coupling to the non- 
strange quarks in the (by a factor of 1.4^) in order to 
reproduce the empirically extracted ii'"''-nucleus interac- 
tion 1^-^,^^, which is slightly repulsive if one wants 
to be consistent with the K^N scattering length, and 
the corresponding value at pB = 0.16 fm^"^ is estimated 
to be about 20 MeV Thus, we will use the stronger 
vector potential, lA'^VJ, for the A'+-meson in this study. 
Calculated mean field potential felt by X^-meson, using 



lA^VJ, is shown in FIG 
usual QMC formalism 
calculate effective masses 




Through Eqs. (|1|) - (|) and 
|7l,|2j ^,^l[| we self-consistently 
m^, and mean field potentials, 
, in symmetric nuclear matter. The scalar (U^) and 

vector (U^) potentials felt by the hadrons, h, in nuclear 
matter are given by: 
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FIG. 1. Total potential JJtot for nucleon and A-hyperon 
shown as a function of the baryon density, ps, in units of 
the nuclear matter saturation density, po=0.15 fm"''. 



In Eq. (|1 1|) , /s is the third component of isospin pro- 
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jection of the hadron, h, and the p meson mean field po- 
tential, , is zero in symmetric nuclear matter. Then, 
within the approximation that the mean field potentials 
are independent of momentum, the four-momentum of 
the hadron is modified by, — i^^/p^ + m]'^ + U^,p), 
which modifies not only the kinematical factors such as 
the flux and the phase space, but also modifies the re- 
action amplitudes. Obviously, the reaction thresholds 
are modified in nuclear matter and now depend on the 
baryon density. 

FIG. 1^ shows the total (Utot) nucleon and A-hyperon 
potentials at zero momenta as function of the baryon den- 
sity, in units of the saturation density of nuclear matter 
po=0.15 fm^'^. Let us recall that at momentum p — 0, 

Utot ^U, + U,, (12) 

and the total potential for the S-hyperon is almost equal 
to that for the A-hyperon. 

FIG. |l| indicates that both nucleon and hyperon po- 
tentials approach minima around normal nuclear matter 
density, which refiects the fact that around po the energy 
density of nuclear matter is minimized. FIG. |^ shows 
the density dependence of the total K and _ftr*-meson 
potentials at zero momenta. The total kaon potential is 
repulsive as explained before, and depends substantially 
on the baryon density. The _ftr*-meson total potential is 
attractive at baryon densities below ~2.7/9o- 
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FIG. 2. Total potential Utot for K+ and if* (892)+ mesons 
plotted as function of the baryon density, ps, in units of sat- 
uration density, po=0.15 fm~^ of the nuclear matter. 



Now, we will discuss the in-medium modification of 
the resonance masses. At present it seems that there 
is no reliable estimate for the in-medium modification of 
masses for the higher mass baryon resonances. In view of 
its numerous successful applications elsewhere, we base 
our estimate on the QMC model |34|j3^ ]. 

We assume that the light quarks in the baryon reso- 
nances are responsible for the mass modification in nu- 
clear medium, as in the QMC model [^,^. However, 
there is a possibility that the excited state light quarks 
may couple differently to the scalar a field from those in 
the ground states, although we expect the difference is 
small. Thus, we estimate the range for the in-medium 
baryon resonance masses by the following two extreme 
cases, i.e., (i) all light quarks including those in the ex- 
cited states play the same role for the mass modifica- 
tion as those in the ground states, (ii) only the ground 
state light quarks play the role as in the usual QMC 
model. These two cases are expected to give the maxi- 
mum and minimum limits for the mass modifications of 
the baryon resonances. Specifically, the range for the 
effective masses of the baryon resonance in medium is 
given (see, e.g., Ref. for the quark model basis of the 
baryon resonances) : 
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These in-medium resonance masses may be expected to 
modify the resonance propagators in the reaction ampli- 
tudes. To avoid introducing extra unknown parameters 
in this initial study, we approximate the in-medium res- 
onance widths appearing in the propagator by the free 
space ones. From Eqs. ( p^ ) - (p7|), we will show results 
for the cross section calculated using the lower limit for 
the resonance masses. However, we have also performed 
the calculation using the upper limit for the resonance 
masses, and confirmed that our conclusion remains the 
same. 



III. RESONANCE MODEL 

We explain in this section the resonance model 
pO|" ^^5H 4^], which could describe the energy depen- 
dence of the total cross sections, ttN YK, quite well, 
and has been used widely in kaon production simulation 
codes [|l^-|l|,||,||,^. We extend the model by includ- 
ing medium modification of the hadron properties, not 
only in the kinematic factors such as the flux and the 
phase space, but also in the reaction amplitudes. 
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We consider kaon and hyperon production processes 
in TriV collisions shown in FIGs. I and |. Because dif- 
ferent intermediate states and final states contribute to 
the TiN AK and TriV — > Eif reactions, the in-medium 
modification of the reaction amplitudes is also expected 
to be different, as will indeed be shown later. 

In TABLE | we summarize the data for the resonances 
which are included in the model. 
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FIG. 3. K and A production processes. 
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FIG. 4. and E production processes. 

The effective Lagrangian densities used for evaluating 
the processes shown in FIGs. and are: 



CttNN = -igTrNNN'-f5fN ■ TT, 



7rJVAf(1650) 



-gTTNN(1650) 



X (TV(1650)fTV • 7? + TVtTV(1650) • tt 



7rJVAr(1710) 



"*.9irAfAf(1710) 



X (TV(1710)75f • 7? -f TV75fA^(1710) • 7? ) 

_ 9iTNN{n20) 
NN{n20) — 



(18) 
(19) 
(20) 



X (Af^(1720)fAr-9^7r + A^fTV''(1720)-9^7f), (21) 



7rAfA(1920) 



57rArA(1920) 



X (a^(1920)2'tV • a^TT + Al''^A^(1920) • .9^7? 

•£iCAAr(1650) = ^5/<AW(1650) 

X (TV(1650)AA: -I- i^AA^(1650)) , 

C-KKNiYlW) = KAN (1710) 

X (TV(1710)75Ais: -f A:A75A^(1710)) , 



C 



KAN(1720) 



9KAN(172a) 

rriK 



X (TV'^(1720)Aa^A: + (a^if)AA-^(1720)) , 

•CKEAf(1710) = ^i9KT.N(1710) 

X (TV(1710)75r • l^K + Rll ■ f75A^(1710) 



C 



KT,N{172Q) 



9KT,N{1720) 



(N''{n20)T-J^df,K + {df,K)J: ■ f TV" (1720)^ 



(22) 

(23) 
(24) 

(25) 
(26) 

(27) 
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A'SA(1920) 



5'KSA(1920) 

rriK 



X (a^(1920)I' • tdf.K + {d^K)l: ■ 2'^A'^(1920)) , (28) 

i^K'-(Sm)KTT = 'i9K'(S92)KTT 

X (^KtK*{892) ■ 9^^- {d''K)TK* (892) ■ cf'j + h.c. 



X {N-f^' AK* (892) + hj 



^K'(892)'SN — ^.9K'(892)S7V 

X (^N-/^'f -l^K* (892) + h 
In the above, the operators 2* and are defined by 



(29) 
(30) 
(31) 



^ (1^-Ml2*^)^ 

£=±1,0 



(32) 
(33) 



£=±1.0 



with M, and Af being the third components of 
the isospin projections, and r the Pauh matrices. 
iV,iV(1710),iV(1720) and A(1920) stand for the fields 
of the nucleoli and the corresponding baryon reso- 
nances. They are expressed by TV = (p, n), sim- 
ilarly for the nucleon resonances, and A(1920) = 
(A(1920)++, A(1920)+, A(1920)°, A(1920)-) in isospin 
space. The physical representations of the kaon field are, 
K'^ = {K+,K^) and K = {K~,K°), respectively, and 
similarly for the ii'*(892), where the superscript T means 
the transposition operation. They are defined as an- 
nihilating (creating) the physical particle (anti-particle) 
states. For the propagators iSp{p) of the spin 1/2 and 
iG^'^{p) of the spin 3/2 resonances we use: 



iSpiv) = i 



'y ■ p + m 



— m? 



im 



Yfuii ' 



(34) 



FK*(S,92)KA\lj^{PK - P^)\) = 

xC|^(w<: -pV)|cxp f -/3|^(pV -P7r)P 



(38) 



In the calculation, the form factors of Eqs. (|37| ) and 
( p8| ) are multiplied by the corresponding coupling con- 
stants. 

TABLE I. Resonances included in the model. Confidence 
levels of the resonances are, A''(1650) * * * *, A''(1710) * **, 
iV(1720) * * * * and A(1920) * ** [Q. Note that the A(1920) 
resonance is treated as an effective resonance which effec- 
tively represents the contributions of six resonances, A (1900), 
A(1905L A(1910), A(1920), A(1930) and A(1940). See 
Refs. |4O|0 for this effective treatment of the A(1920). 



Resonance 


Width Channel B.R. ratio Used 


N(1650) (i") 


150 N-K 0.60 - 0.80 0.700 
(MeV) KK 0.03 - 0.11 0.070 


N(1710) (i^) 


100 Ntv 0.10 - 0.20 0.150 
(MeV) KK 0.05 - 0.25 0.150 
T^K 0.02 - 0.10 0.060 


N(1720)(r) 


150 Ntv 0.10 - 0.20 0.150 
(MeV) KK 0.03 - 0.10 0.065 
Sif 0.02 - 0.05 0.035 


A(1920) (1+) 


200 Ntv 0.05 - 0.20 0.125 
(MeV) T,K 0.01 - 0.03 0.020 


K*{m2) (1-) 


50 Ktv ~ 1.00 1.00 

(MeV) 



ic^ip) = i- 



-pf"'{p) 



(35) 



with 

P^^^ip) 



{-f ■ p + m) 
1 



77^7" 



3 m 



3m2 



(36) 



where m and F-''"" stand for the mass and full decay 
width of the corresponding resonances. For the form fac- 
tors, F{q) (q is the momentum of meson, tt or K), ap- 
pearing in the meson-baryon- (baryon resonance) vertices, 
we use: 



m = 



A2 



(37) 



For the K* {892)-K-7t vertex we adopt the form factor of 
Rcf. M: 
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The parameters of the model, namely the form factors 
in the interaction vertices and the coupling constants, 
were fixed by available experimental data on the different 
ir+N^Y+K reaction channels. Furthermore, the same 
parameters which were determined by the n+N—^Y+K 
reactions were used in the calculations of strangeness pro- 
duction in baryon-baryon collisions, and t hey also repro- 
duced the available data reasonably well |^5|-^ . In TA- 
BLE |l| we list all values for the coupling constants, cut- 
offs, C and (3 parameters relevant for the present study. 

TABLE IL Values for coupling constants, cut-offs, C and 
l3 parameters used in the present study. For details about 
the coupling constants relevant for A(1920), g^NA(i920) ^nd 
ffifEA(i920)) see Refs. 



vertex 




cut-off A (MeV) 


TriViV 


UA 


1050 


nNN{l650) 


1.12 X 10"^ 


800 


nNNlnW) 


2.05 X 10"^ 


800 


7riVA^(1720) 


4.13 X 10"^ 


800 


7rArA(1920) 


1.13 X 10-1 


500 


Ji:AiV(1650) 


5.10 X 10"^ 


800 


KAN{niO) 


3.78 


800 


KAN{1720) 


3.12 X 10-1 


800 


7<'E7V(1710) 


4.66 


800 


A'EiV(1720) 


2.99 X 10"! 


800 


Js'EA(1920) 


3.08 X 10-1 


500 


J^*(892)AiV 


1.62 X 10-^ 


1200 


K*{892)T.N 


1.62 X 10-^ 


1200 


K*{892)Kn 


5.48 X 10-^ 


C = 2.72fm 
P = 8.88 X IQ-'^fm^ 



Since the strength of the KK and Y.K coupling to the 
various baryonic resonances are different (see TABLE 1), 
the dynamics of A and S production are also different. 
This can be understood as follows. 

The amplitude of the resonance propagator in the re- 
action amplitude, ■nN — )■ Y K, becomes maximal when 
the invariant collision energy crosses the mass of the res- 
onance. Furthermore, the cross section a{nN — s- YK) 
is proportional to (qy)^^'^^ near threshold (with £ and 
qy — IoyI the orbital angular momentum and momen- 
tum of the YK pair in the center of mass system). In 
seeking to understand the difference in the behaviour of 
a{7TN KK) and a{'KN TiK) it is important to note 
that the former receives a large contribution from the s- 
wave Af(1650) (with £ = 0), while the latter is dominated 
by p-wave resonances (with ^ = 1). The subtle interplay 
between the change in momentum dependence associated 
with threshold variations and the change in resonance 
amplitudes associated with varying resonance masses is 
responsible for the different behaviour of (t{ttN KK) 
and cr(7r7V T,K). 

The total free mass of the YiK system, mj] -I- rrij^, is 
very close to the mass of the Pii(1710) resonance, i.e., 
Mpjj — ms — rriK —27 MeV, which dictates its substan- 
tial role in E hyperon production The situation is 
different for the KK system, since the nearest baryonic 
resonance is S'ii(1650) and Ms^^ — ttia — rriK —41 MeV. 
Furthermore, because S'ii(1650) and Pii(1710) have dif- 
ferent couplings to the final states, KK (relative s-wave) 
and Y.K (relative p-wave) , respectively, the change in the 
vertices, S"!! (1650) AX and Pii(1710)SA' due to the mo- 
mentum modified in medium are also different. Further- 
more, when threshold changes as the baryon desnsity 
varies, the change in the strength of the contributions 
from the 511(1650) to the ttN KK reaction and that 
of the P33(19 20) to the t:N SAT reaction should be 
different as a function of baryonic density. Thus, com- 
bining these effects for the ttN KK and ttN — > EAT 
reactions, we can see that the energy dependence of the 
cross sections in a nuclear medium can be significantly 
modified. 



IV. THE n+N^Y+K+ REACTION THRESHOLD 
IN NUCLEAR MATTER 

The dispersion relation in nuclear matter relating the 
total energy E and the momentum p of the particle is 
written as 

E=^p-' + {m + Us)^ + U,, (39) 

where m is the bare mass and Us, Uy denote the scalar 
and vector parts of the potential in nuclear matter. The 
threshold, y^sth, for the reaction n+N^Y+K'^ is given 
as the sum of the total energies of the final Ar+-meson 
and y-hyperon, taking their momenta to be zero and 
hence 
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yi^ = m^ + m^ + [/f +[// + C/J, (40) 

where now the upper indices denote kaons and hypcr- 
ons. The soHd Hnes in FIG. ^ show the K and X+S 
reaction thresholds, ^/sth, as a function of the baryon 
density. Obviously, in free space the scalar and vector 
potential vanish and the reaction threshold equals to the 
sum of the bare masses of the produced particles, which 
is shown by the dashed lines in FIG. js] for the A. and 
K^Yi final states. 




Pb/Po 

FIG. 5. The threshold energy, ^/s^, for K+ A and if+E 
production given by their total in-medium energy at zero mo- 
mentum, as a function of the baryon density, ps, in units 
of saturation density of nuclear matter, po=0.15 fm~^. The 
solid lines indicate our results, while the dashed lines show 
the threshold in free space. 



also found that at baryon densities pB>0.2 fm~'^ the 
strangeness production threshold in TriV collisions is 
higher than the free space case. 

V. THE TT+N^A+K REACTION IN NUCLEAR 
MATTER 

Now we apply the resonance model to calculate the 
in-medium amplitudes. We keep the coupling constant 
as well as the form factors at the values found in free 
space. While this assumption certainly cannot be com- 
pletely correct in nuclear matter, there are no presently 
established ways to improve this part of our calculation. 
In principle, since we started from the reaction amplitude 
itself, it is possible to include in-medium modifications of 
the coupling constants as well as the form factors when 
reliable calculations of the changes of these quantities in 
nuclear matter become available. In the following calcu- 
lations we include the vector and scalar potentials for the 
interacting (initial) nucleons and final kaons and hyper- 
ons, as well as for the intermediate baryonic resonances 
and i^*-meson. 

FIG. ^ shows our results for the differential cross sec- 
tion for the reaction 7r^+j3^A+i^" at the invariant col- 
lision energy y^=1683 MeV. It is calculated both in free 
space (the solid line) and in nuclear matter, at baryon 
density ps—Po (the dashed line) and pB='^po (the dot- 
ted line). For comparison, we also show in FIG. ^ the 
experimental data collected in free space |^,|5^. The 
important finding is that not only the absolute magni- 
tude, but also the shape (the dependence on the cos6') of 
the 7r~ -l-p^ A-|-i^° differential cross section, depends on 
the baryon density. 



It is important, that while the /C^-meson energy at 
zero momentum increases with the baryon density (see 
FIG. ^), because of the negative A and S potentials 
the reaction threshold in nuclear matter at pB<1.4po is 
shifted below that in free space. 

The maximal downward shift of the reaction thresh- 
old in nuclear matter occurs at baryon densities around 
PB~0.6po. This value is the result of competition be- 
tween the simple, linear dependence on density of the 
vector potentials and the more complicated, non-linear 
behaviour of the scalar potentials. (A similar competition 
leads to the saturation of the binding energy of normal 
nuclear matter in the QMC model.) Furthermore, the 
maximum of the downward shift of the ■k-\-N^Y+K~^ 
reaction threshold amounts to roughly 30 MeV. We 
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FIG. 6. The yr'+p^A+AT^ differential cross section at in- 
variant collision energy y^=1683 MeV as a function of the 
cosine of the kaon emission angle in the center of mass sys- 
tem. The experimental data are from Ref. (the squares) 
and from Ref. jsoj (the circles). The lines show our calcu- 
lations in free space (solid) and in nuclear matter at baryon 
density pB=po (dashed) and pB—2po (dotted), with po = 0.15 
fm-^ 



One of the simplest ways to construct the in-medium 
reaction cross section is to take into account only the in- 
medium modification of the flux and phase space factors 
while leaving amplitude in matter the same as that in 
free space, without including medium effect To shed 
more light on the problem of how the reaction amplitude 
itself is modified in nuclear matter, we show in FIG. |^ 
the reaction amplitudes squared in arbitrary units for 
the TT^+p^A+K^ reaction, calculated at yG—lJ GeV 
and 1.9 GeV, in free space (the solid lines), pb — Po 
(the dashed lines) and ps = 2po (the dotted lines). Our 
calculation clearly indicates that the tt^ +p^A+K'^ re- 
action amplitude in nuclear matter differs substantially 
from that in free space at these energies, and that the 
amplitudes depend on the baryon density. 
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FIG. 7. The (dimensionless) invariant amplitude squared 
for the tt' +p^A+K^ reaction, as a function of cos 6 (the 
isr+-meson emission angle in the center of mass system), cal- 
culated for the invariant collision energies -^8=1.7 GeV (up- 
per) and 1.9 GeV (lower). The lines show the result for free 
space (solid) and for nuclear matter at baryon density pB~po 
(dashed) and pB~2po (dotted). 



Finally, 

the energy dependence of the total tt^ +p^A+K'^ cross 
section is shown in FIG. |^, as a function of the invariant 
collision energy, ^/s. The experimental data in free space 
are taken from Ref. The calculations for free space 
are in reasonable agreement with the data, as shown by 
the solid line. The dashed line in FIG. H shows the results 
obtained for nuclear matter at pB=Po, while the dotted 
line is the calculation at pB=2po- 

Clearly the total n^'+p^A+K'^ reaction cross section 
depends substantially on the baryon density. Further- 
more, as already discussed in Sec. Ill, the reaction thresh- 
old at baryon density pB=po is shifted downward as com- 
pared to that in free space, while at pB='2po it is shifted 
upwards. 
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FIG. 8. Energy dependence of the total cross section, 
TT^+p^A+Ji''', as a function of the invariant colhsion energy, 
^/s, calculated for different baryon densities. The data in free 
space are taken from Ref. [^. The lines indicate our results 
for free space (solid) and for nuclear matter at baryon density 
Pb=Po (dashed) and pB=2po (dotted). (Only the solid curve 
should be compared directly with the data.) 



Obviously, heavy ion collisions probe a range of baryon 
densities from pb=0 up to several times normal nu- 
clear matter density, po- The calculation of the time 
and spatial dependence of the baryon density distribu- 
tion is a vital aspect of dynamical heavy ion simula- 
tions. However, a first estimate of the density averaged 
total 7r~+p^A+K^ cross section can be gained from 
FIG. ||. Of course, the data is only available in free 
space and should only be directly compared with the solid 
curve. Nevertheless, it is suggestive for the problem of in- 
medium production to note that a crude average of the in- 
medium cross sections over the range 0<pB<2po would 
be quite close to the free space cross section at energies 
around the free space threshold. This seems to provide a 
reasonable explanation of why the heavy ion calculations 
including [|l^,^^ the n+p^A+K cross section in free 
space, that is without a repulsive kaon potential, can re- 
produce the data |ll6|-po|. More quantitative calculation 
and discussion of this effect will be given in Sec. VI. 



VI. THE n+N- 



>T.+K REACTION IN NUCLEAR 
MATTER 



The TT+N-^Yi+K reaction involves different dynamics 
in comparison with the n+p^A+K reaction, because the 



reaction involves the different intermediate baryonic res- 
onances. For instance, although the iV(1650) resonance 
couples to AiV channel, it does not couple to the 'SN 
state. Moreover, the A(1920) resonance couples to UN 
channel, but does not couple to the AA^ channel in our 
model. For this reason, the dependence on the baryon 
density of the reaction in nuclear matter, TT+N—fTi+K , 
is quite different from that of the n+N^A+K. 




FIG. 9. Energy dependence of the total cross section, 
tt'^ +p^'E'^ , as a function of the invariant collision en- 
ergy, -^s, calculated for different baryon densities. The data 
in free space are taken from Ref. |Q. The lines indicate our 
calculations for free space (solid) and for nuclear matter at 
baryon density pB=po (dashed) and pB~2po (dotted). 



In FIG. I we show the energy dependence of the to- 
tal cross section, 7r++p^E+-|-if +, as a function of the 
invariant collision energy, y^. The experimental data in 
free space are taken from Ref. [|5^. The free space data 
are well reproduced by the calculations in free space, as 
shown in FIG. ^ by the solid line. The dashed line indi- 
cates the results obtained for nuclear matter at baryon 
density pB=Po, while the dotted line shows the result at 
PB='2po- 

Again, as already discussed in Sec. Ill, the density de- 
pendence of the hadron masses and the vector potentials 
leads to a shift of the reaction thresholds in nuclear mat- 
ter. Because of the density dependence of the E-hyperon 
potential, the threshold at normal nuclear matter den- 
sity {pb=Po) is shifted downwards compared with that 
in free space. At pB—'2po the 'SK reaction threshold is 
shifted upwards relative to the threshold in free space. 
Moreover, the magnitude of the tt~^+p^Y,^+K~^ cross 
section depends much more strongly on the density than 
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the 7r~+p^A+ivr" reaction. 

FIGs. |l^ and [l^ show the energy dependence of 
the total cross sections for the 7r~+p-^I]"+i^° and 
TT" +p^Y,~ +K'^ reactions, respectively. The data in free 



space 1 52 are well reproduced with the calculations in 
free space, which are indicated by the solid lines. The 
cross sections calculated for the nuclear matter, except 
for 7r~+p^I]°+_ftr" at pB = 2po, are substantially en- 
hanced in comparison with those in free space, at energies 
above near the in-medium reaction thresholds. 




FIG. 10. Energy dependence of the total cross section, 
Tv~ +p^TP +K'' , as a function of the invariant collision en- 
ergy; y/s, calculated for different baryon densities. The data 
in free space are taken from Ref. The lines indicate our 
calculations for free space (solid) and for nuclear matter at 
baryon density pB~po (dashed) and ps=2po (dotted). 




FIG. 11. Energy dependence of the total cross section, 
7r~4-p— >E~+A'+, as a function of the invariant collision en- 
ergy, ^/s, calculated for different baryon densities. The data 
in free space are taken from Ref. The lines indicate our 
calculations for free space (solid) and for nuclear matter at 
baryon density pB=po (dashed) and pB—^po (dotted). 



VII. IMPACT ON HEAVY ION STUDIES 

It is expected that in relativistic heavy ion collisions 
at SIS energies nuclear matter can be compressed up to 
baryonic densities of order /9B~3po The baryon den- 
sity pb available in heavy ion collisions evolves with the 
interaction time, t, and is given by the dynamics of the 
heavy ion collision. Moreover, the density is large in the 
very center of the collision. In the following estimates 
we investigate the density dependence of the production 
cross section for central central heavy ion collisions. How- 
ever, it should be remembered that at the edges, where 
most particles are expected to be located, the density 
dependence of the strangeness production mechanism is 
not strong compared to that of the central zone of the 
collision. 

To calculate the 7^+-meson production cross section 
averaged over the available density distribution we adopt 
the density profile function obtained by dynamical sim- 
ulations for Au+Au collisions at 2 AGeV and at 
impact parameter b=0. This can be parametrized as 



PB{t) = Pn 



exp 



(41) 



where the parameters, Pmax='^Poi t=13 fm and 
At=6.7 fm, were fitted to the heavy ion calculations [p3[. 
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is given in Ref. |Q as a function of the invariant collision 
energy, ^/s. It can be parametrized for GeV as 




FIG. 12. Energy dependence of the total cross section, 
TT^+p^A+A'", as a function of the invariant collision energy, 
y/s. The data in free space are taken from Ref. |Q . The solid 
line indicates our calculation for free space. The dashed line 
shows the cross section calculated by averaging over the den- 
sity function profile jHsj given by the time evolution obtained 
for Au+Au collisions at 2 AGeV (see Eq. (|l|)). 



The total cross section for the +p^K+K^ reaction 
integrated over the time range 5<t<23 fm and weighted 
by the time dependent density profile given in Eq. (41), 
is shown by the dashed line in FIG. |lj. The limits of 
the t integration were taken from the simulations of the 
Au+Au collision time evolution in Ref. (ss). The circles 
and solid line in FIG. |l^ show the experimental data in 
free space and the calculations in free space, respec- 
tively. 

One can see that the total cross section averaged over 
the collision time (time dependent density profile) for the 
TT~+p—^A+K° reaction is quite close to the result given 
in free space integrated up to at energies above the pro- 
duction threshold, up to y/s ~ 1.7 GeV. That the results 
shown in FIG. ^ actually explain why the heavy ion cal- 
culations with the free space kaon production cross sec- 
tion might quite reasonably reproduce the experimental 
data, will be discussed more quantitatively below. 

As a matter of fact, the total cross section averaged 
over the time dependent density profile, shown by the 
dashed line in FIG. |l^, should additionally be averaged 
over the invariant collision energy distribution available 
in heavy ion reactions. The number of meson-baryon col- 
lisions, N„iB, for the central Au+Au collisions at 2 AGeV 



dN, 



mB 



dy/s 



No exp 



(42) 



where the normalization factor A'o=6xl0* GeV ^, while 
GeV and A^=0.63 GeV. Note that, at SIS en- 
ergies N„iB is almost entirely given by the pion-nucleon 
interactions, and heavy meson and baryon collisions con- 
tribute only to the high energy tail of the distribution 
in Eq. - with quite small densities |Q. Finally, if 
we also average the calculated, in-medium, total cross 
section for ir^ +p^A+K'^ , shown by the dashed line in 
FIG. O, over the available energy distribution given in 
Eq. (^2^, we obtain an average total kaon production 
cross section of <K>=65 /ib for central Au+Au colli- 
sions at 2 AGeV. This result is indeed compatible with 
the calculations using the free space total cross section 
of the Tr~ +p—>A+K^ reaction, which provide an aver- 
age total kaon production cross section of <K>~71 fjh 
for central Au+Au collisions at 2 AGeV. Note that the 
inclusion of even a slight modification of the K'^ mass be- 
cause of the nuclear medium (without the corresponding 
changes introduced here) leads to a substantial reduction 
of the inclusive spectra (by as much as a factor of 2 
or 3), compared to that calculated using the free space 
properties for the relevant hadrons 1 14 1 . 

We stress that at SIS energies reaction channels with a 
S-hyperon in the final state play a minor role, because of 
the upper limit of the energy available in the collisions. 
As was illustrated by FIG. |^, the downwardly shifted 
TT+N^J^+K reaction threshold at small densities is still 
quite high compared to that for the reaction with the 
A-hyperon in the final state. 



VIII. SUMMARY 

We have calculated the in-medium modification 
of kaon production in pion-nucleon collisions in nu- 
clear matter using the resonance model developed in 
Refs. po|-^. To evaluate the in-medium ii'-meson pro- 
duction for the reaction channels with A and S hy- 
perons in the final states, the density dependent po- 
tentials for the initial, final and intermediate mesons 
and baryons (resonances) were introduced to the reso- 
nance model amplitudes. The vector and scalar poten- 
tials were calculated within the quark-meson coupling 
model @J2|,|^-||. The tt+A^^A+A: and vr+A^E+A 
cross sections were calculated for different baryon densi- 
ties of the nuclear matter. We found, that not only are 
the initial flux and the final phase space of the reactions 
modified in baryonic matter, but the reaction amplitudes 
themselves are also modified. 

It was shown, that the total tt+N^A+K and 
TT+N^T,+K cross sections depend substantially on the 
baryonic density. Furthermore, the reaction thresholds 
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and the absolute magnitudes as well as the dependence 
of the production cross section on the invariant collision 
energy, ^/s, all vary strongly with the density of the nu- 
clear matter. 

To evaluate the impact of our microscopic calculations 
on the heavy ion results, we averaged the kaon produc- 
tion cross section over the baryon density profile, which 
depends on the evolution time of the heavy ion collision. 
Furthermore, in order to compare with the experimental 
data more quantitatively, we calculated the effective to- 
tal kaon production cross section by averaging over the 
invariant collision energy distributions available in heavy 
ion reactions. We found that at low collision energies, 
the density or time averaged ii'^-meson production total 
cross section, calculated using the in-medium properties 
for the meson, hyperons and relevant hadrons, is very 
close to that calculated using the total cross section given 
in free space. 

Thus, our results provide an explanation of why the 
analyses ||l^,|l5|,|l 22 1 of available data on produc- 
tion from heavy ion collisions at SIS energies 18-2l|] 
for the K'^ spectra, can be reasonably described when 
one neglects any in-medium modification of the kaon and 
hadronic properties - i.e., adopting the if "''-meson pro- 
duction cross section given in free space. 

In conclusion, our present study shows that if one ac- 
counts for the in-medium modification of the production 
amplitude (i.e., the in- medium properties of the K'^- 
meson and hadrons) correctly, it is possible to under- 
stand production data in heavy ion collisions at SIS 
energies, even if the if^-meson feels the theoretically ex- 
pected, repulsive mean field potential. The apparent fail- 
ure to explain the production data if one includes the 
purely kinematic effects of the in-medium modification 
of the iir"'"-meson and hadrons, appears to be a conse- 
quence of the omission of these effects on the reaction 
amplitudes. 
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